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SU^WARY 


Estimates  of  the  test  section  wall  interference  effects  have  been 
computed  for  an  elliptic  missile  body  configuration  in  the  Aerodynamic 
Wind  Tunnel  (4T)  of  the  Propulsion  Wind  Tunnel  Facility.  Results  are 
presented  for  test  Mach  numbers  of  0.95  and  1.05.  Comparisons  of  test 
section  wall  static  pressure  measurements  in  4T  with  computations  show 
relatively  good  agreement  at  Mach  0.95  and  follow  the  general  trend  at 
Mach  1.05.  The  wall  interference  produces  a  2.5%  decrease  in  lift 
coefficient  at  the  highest  angle  of  attack  for  which  calculations  were 
made. 
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NOMENCLATURE 


Cl 

Cm 


M 


Lift  coefficient 

Pitching  moment  coefficient 

Pressure  coefficient  (p-Poo)/qoo 

Interference  pressure  coefficient  (Cp  unconfined  flow 
-  Cp  confined  flow) 

Mach  number 


P 

q 

a 

ACl 


ACm 

ACp 

Subscripts 


Pressure,  psfa 

Dynamic  pressure,  psf 

Model  angle  of  attack,  degrees 

Interference  lift  coefficient  (Cl  unconfined  flow 
-  Cl  confined  flow) 

Interference  pitching  moment  coefficient 
(Cm  unconfined  flow  -  C^  confined  flow) 

Wall  static  pipe  Cp  at  lift  coefficient  of 
interest  -  wall  static  pipe  Cp  at  zero  lift 


Test  section  reference  conditions 
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1.0  INTRODUCTION 


Wall  Interference  produced  by  a  tunnel  upon  a  model  is  a  contri¬ 
butor  to  errors  inherent  to  measured  wind  tunnel  data.  In  every  tran¬ 
sonic  wind  tunnel  test,  wall  interference  exists  to  some  degree.  There¬ 
fore,  it  is  advantageous  to  correctly  identify  the  level  of  wall  inter¬ 
ference  in  wind  tunnel  tests. 

Recently  an  elliptic  missile  configuration  (Fig.  1)  was  tested  in 
Tunnel  4T  to  provide  a  data  base  for  verification  of  computational  fluid 
dynamics  (CFD)  methods.  Because  of  the  relatively  large  model  size 
{1.2%  blockage)  and  high  transonic  Mach  numbers  tested,  wall 
interference  effects  were  of  major  concern.  The  ARO-1  code,  an  Euler 
solver,  developed  at  AEDC  (Ref.  1)  is  used  in  this  analysis  to  estimate 
the  wall  interference  on  the  elliptic  missile  configuration.  The  method 
consists  of  computing  the  flow  field  about  a  body  within  a  tunnel  with 
the  porous  walls  modelled  by  an  empirical  correlation  and  then  computing 
the  unconfined  flow  about  the  same  body.  The  wall  interference  is  the 
difference  between  the  two  solutions. 

The  calculations  were  made  on  the  elliptic  body  at  free-stream  Mach 
numbers  of  0.95  and  1.05.  The  angle  of  attack  was  10  deg  for  both  Mach 
numbers;  a  flow  case  of  a  =  20  deg  and  Moo  =  0.95  was  also  computed.' 

Calculations  of  wall  interference  effects  on  an  elliptic  missile 
body  geometry.  Fig.  1,  was  requested  by  Air  Force  Wright  Aeronautical 
Laboratories  (AFWAL/FIMG) ,  Wright-Patterson  AFB,  Ohio,  under  Program 
Element  62201F,  Control  Number  2404.  The  program  manager  was  Lt.  Norman 
Schmoeker,  AEDC/DOFA.  The  results  were  obtained  by  Calspan 
Corporation/AEDC  Division,  operating  contractor  for  the  aerospace  flight 
dynamics  testing  effort  at  the  AEDC,  AFSC,  Arnold  Air  Force  Station, 
Tennessee.  The  analysis  was  conducted  in  the  Propulsion  Wind  Tunnel 
Facility  (PWT)  during  the  period  from  17  September  through  20  December 
1985  under  AEDC  project  number  CD48PB,  PWT  Test  Number  TC793.  Docu¬ 
mentation  of  the  test  is  contained  in  Ref.  2. 

2.0  COMPUTATIONAL  MESH 

The  process  of  computing  wall  interference  in  the  present  method¬ 
ology  requires  two  computational  meshes.  One  mesh  models  the  space 
between  the  aerodynamic  body  and  the  tunnel  walls  to  calculate  the  con¬ 
fined  flow  in  the  tunnel.  This  simulation  is  done  by  applying  approp¬ 
riate  porous  wall  boundary  conditions  at  the  mesh  outer  boundary.  Com¬ 
putational  space  is  then  added  to  the  exterior  of  the  first  mesh  resul¬ 
ting  in  the  second  mesh  for  calculations  of  the  unconfined  flow  about 
the  body.  With  the  additional  space  and  freestream  conditions  imposed 
on  the  second  mesh  outer  boundary,  unconfined  flow  on  the  body  may  be 
simulated.  The  difference  of  the  two  solutions  is  the  wall  inter¬ 
ference.  Two  meshes  are  required  for  each  angle  of  attack  case  con¬ 
sidered  because  of  the  change  in  model  attitude  relative  to  the  tunnel 
walls. 
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The  B25  elliptic  body  configuration  (Fig.  1)  was  described  ana¬ 
lytically  and  a  computational  surface  grid  constructed.  The  present 
method  requires  the  body  to  be  closed;  therefore,  a  conical  section  was 
added  to  the  base  of  the  body.  This  addition  simulates  the  wake  region 
better  than  a  flat  base.  Because  of  symmetry,  only  half  of  the  tunnel 
and  model  were  modelled  to  minimize  the  mesh  points  required. 
Approximately  900  mesh  points  were  used  on  half  the  model  surface. 

The  surface  grid  was  transformed  into  a  warped  spherical  grid 
topology  to  be  used  as  the  inner  boundary  of  the  tunnel  mesh.  Two  views 
of  the  surface  mesh  are  shown  in  Fig.  2.  After  analytically  describing 
the  tunnel  surface,  the  inner  model  surface  and  outer  tunnel  surface 
were  used  as  boundary  conditions  for  a  three-dimensional  Poisson 
equation  solver  (Ref.  3)  to  create  the  interior  mesh.  Two  planar  views 
of  the  model  at  10°  angle  of  attack  are  shown  within  Tunnel  4T  (Figs.  3 
and  4)  to  demonstrate  the  interior  mesh.  Also,  the  model  surface  is 
shown  positioned  inside  the  tunnel  surface  in  Fig.  5. 

3.0  RESULTS  AND  DISCUSSION 

Computations  of  wall  interference  have  been  made  on  the  B25 
elliptic  missile  body  configuration  in  Tunnel  4T.  The  calculations  were 
obtained  at  three  flow  conditions: 

1)  Moo  =  0.95,  a  =  10  deg 

2)  Moo  =  1.05,  a  =  10  deg 

3)  Moo  =  0.95,  a  =  20  deg 

The  wall  porosity  for  all  cases  is  5%.  The  results  are  presented  sepa¬ 
rately  for  each  of  the  above  cases.  In  each  case,  the  wall  interference 
is  illustrated  as  a  pressure  coefficient  increment  produced  by  the  wall. 
Therefore,  the  pressure  contours  shown  in  Figs.  6-14  are  the  differences 
(ICp)  between  the  unconfined  flow  pressures  and  the  tunnel  confined  flow 
pressures.  In  addition,  the  pressure  differences  were  integrated  over 
the  body  to  obtain  the  incremental  interference  lift  and  pitching 
moment  coefficients.  These  coefficients  are  presented  in  Table  1  for 
the  three  cases  investigated. 

As  illustrated  in  Fig.  15,  the  ARO-1  code  underpredicted  the  lift 
on  the  model;  consistently  lower  lift  curve  slopes  are  always  evident  in 
the  calculations  when  compared  to  experimental  data.  The  measured  lift 
coefficient  data  were  obtained  from  the  Elliptic  Missile  Body  Force  Test 
documented  in  Ref.  4.  The  mismatch  occurs  because  viscous  effects  are 
ignored  and  computer  storage  restrictions  do  not  allow  complete 
resolution  of  the  entire  model.  However,  the  results  are  usable  if  one 
recognizes  that  the  increment  attributable  to  wall  interference  is 
proportional  to  lift.  Therefore,  the  comparisons  with  experimental  data 
should  be  made  with  matched  lift  coefficients  rather  than  angle  of 
attack. 
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3.1  Moo  =  0.95,  a  =  10  deg 


The  pressure  difference  distribution  on  the  elliptic  body  is  shown 
in  Fig.  6  for  the  top  and  bottom  surfaces.  The  wall  interference  on  the 
body  is  initially  small  near  the  nose  but  increases  axially  along  the 
body.  Also,  the  bottom  surface  pressure  difference  contours  are  shifted 
aft  from  the  corresponding  contours  on  the  top  surface  illustrating  the 
angle  of  attack. 

The  wall  interference  effect  at  the  tunnel  walls  is  shown  in  Fig. 
7.  The  proximity  of  the  base -of  the  body  to  the  bottom  wall  is  evident 
by  the  high  gradient  and  level  of  contours  on  the  wall.  At  the  tunnel 
symmetry  plane,  the  patterns  of  pressure  difference  in  the  tunnel  are 
shown  in  Fig.  8.  The  pressure  difference  contours  on  the  base  conical 
addition  are  significant,  but  since  this  addition  does  not  exist  in  the 
actual  testing  of  the  model,  these  contours  may  be  disregarded. 


3.2  Moo  =  1.05,  a  =  10  deg 

The  effect  of  a  higher  free-stream  Mach  number  is  shown  in  this 
flow  case.  The  pressure  difference  contours  are  at  a  lower  absolute 
level  on  the  body  (Fig.  9)  however,  the  interference  lift  coefficient  is 
higher  as  shown  in  Table  1.  The  pressure  difference  contours  are  at  a 
higher  absolute  level  on  the  tunnel  walls  (Fig.  10)  with  the  same  basic 
patterns  as  Case  1.  The  interference  effects  at  the  symmetry  plane 
(Fig.  11)  are  also  consistent  with  the  previous  case  but  show  a  more 
pronounced  interference  pattern  in  the  tunnel. 

3.3  Moo  =  0.95,  a  =  20  deg 

This  flow  case  shows  the  effect  of  a  higher  angle  of  attack  on  the 
wall  interference  in  the  tunnel.  As  shown  in  Figs.  12  thru  14,  the 
pressure  difference  contour  levels  have  intensified  on  the  body  and  the 
tunnel  walls.  The  contour  patterns  have  also  changed  on  the  tunnel 
walls  since  the  model  attitude  is  different.  The  interference  is  exten¬ 
sive  as  compared  to  the  previous  cases. 

Considering  these  three  flow  cases,  the  results  show  that  a  higher 
free-stream  Mach  number  causes  higher  levels  of  wall  interference  on  the 
model  based  on  interference  lift  coefficient.  Also,  higher  angles  of 
attack  result  in  a  significant  increase  in  wall  interference. 

4.0  WALL  PRESSURE  COMPARISONS 

In  this  section,  wall  pressure  measurements,  taken  by  static 
pressure  pipes  located  on  the  centerline  of  the  top  and  bottom  walls, 
are  compared  to  calculations  on  the  top  and  bottom  walls  for  four  cases: 

1)  Moo  =  0.95,  Cl  =  0.6 

2)  Moo  =  0.95,  Cl  =  2.1 

3)  Mo,  =  1.05,  Cl  -  0.6 

4)  Mo,  =  1.05,  Cl  =  2.1 
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Because  of  local  effects  of  the  static  pipe,  wall  pressure  data  at  zero 
lift  were  subtracted  from  wall  pressure  data  of  each  case  considered  and 
compared  to  calculations  made  similarly.  These  pressure  differences 
(ACp)  are  plotted  versus  tunnel  station  in  Fig.  16. 

The  wall  pressure  comparisons  at  Moo  =  0.95  and  Ci_  =  0.6  (Fig.  16a) 
show  good  agreement  on  the  top  wall  with  fair  prediction  of  the  pressure 
trends  on  the  bottom  wall.  For  the  higher  lift  case  (Fig.  16b),  the 
same  trends  occur  on  both  walls.  For  the  Moo  =  1.05  cases  (Figs.  16c  and 
16d),  the  data  is  scattered  enough  that  no  good  comparison  can  be  made. 
However,  the  calculations  show  a  consistent  trend  as  compared  to  the  Mo, 
=  0.95  cases  and  a  general  trend  may  be  seen  in  the  data  as  compared  to 
the  calculations.  This  data  scatter  may  occur  because  of  inherent 
problems  associated  with  testing  near  Moo  =  1. 

Based  on  previous  experience,  the  quality  of  the  comparisons 
between  calculated  and  measured  pressure  coefficients  near  the  wall  is 
an  indication  that  the  actual  wall  effects  on  the  model  should  be  rela¬ 
tively  close  to  those  of  the  calculations. 

5.0  CONCLUSIONS 

Computations  of  wall  interference  on  the  B25  elliptic  missile  con¬ 
figuration  have  been  presented.  Agreement  between  the  measured  and 
computed  wall  pressures  provides  a  level  of  confidence  that  the  computed 
interference  on  the  body  is  good.  The  wall  interference  produces  a  2.5% 
decrease  in  lift  coefficient  at  the  highest  angle  of  attack  for  which 
calculations  were  made.  These  wall  interference  effects  should  be 
accounted  for  if  the  measured  data  are  to  be  used  to  verify 
computational  fluid  dynamics  codes. 
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Figure 


a.  Top  View 

Computational  Model  Surface 


b.  Side  View 


Figure  2.  Concluded 


Figure  3 


Computational  Mesh  in  Horizontal  Plane 


Figure  4,  Computational  Mesh  in  Symmetry  Plane 
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b.  Bottom  Surface 
Figure  6.  Concluded. 
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Figure  7.  Concluded. 
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Figure  8.  Symmetry  Plane  Contours  of  Pressure  Differences, 

M  =  0.95,  a  =  10  deg. 
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a.  Top  Surface 

Figure  9.  Distribution  of  Interference  on  Model  Surface, 

M  =  1.05,  a  =  10  deg. 
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b.  Bottom  Surface 
Figure  9.  Concluded. 
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a.  Top  and  Side  Surface 

Figure  10.  Distribution  of  Interference  on  Tunnel  Boundary, 

=  1.05,  a  =  10  deg. 
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b.  Bottom  and  Side  Surface 
Figure  10.  Concluded. 
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Figure  11.  Symmetry  Plane  Contours  of  Pressure  Differences, 

M  =  1.05,  a  =  10  deg. 
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a.  Top  Surface 

Figure  12.  Distribution  of  Interference  on  Model  Surface, 
M^  =  0.95,  a  =  20  deg. 


b.  Bottom  Surface 
Figure  12.  Concluded. 
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a.  Top  and  Side  Surface 

Figure  13.  Distribution  of  Interference  on  Tunnel  Boundary, 

=  0.95,  a  =  20  deg 
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b.  Bottom  and  Side  Surface 
Figure  13.  Concluded. 
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Figure  14. 


Symmetry  Plane  Contours  of  Pressure  Differences, 
M  =  0.95,  a  =  20  deg. 
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Figure  15,  Computed  and  Measured  Lift  vs  Angle-of-Attack 
in  PWT-4T,  M  «  0.95 
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a.  M  =  0.95,  Cl  =  0,6 

Figure  16.  Computed  and  Measured  Wall 
Comparisons  in  PWT-4T 
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b.  M  =»  0,95,  Cl  »  2.1 
Figure  16.  Continued 
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Figure  16.  Continued 
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d.  M  =  1.05,  Cl  =  2.1 
Figure  16.  Concluded 
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Table  1.  Summary  of  Integrated  Wall  Interference  Effects 
on  Model  Force  and  Moment 


Cl  =  0 

Cpi  =  0 

ACl 

0.0177 

AC[p 

0.0033 

Moo  =  0.95 

6 
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Cl  =  2.1 
=  0.86 
0.0540 
-0.0002 


Moo  =  1.05 


Cl  =  0 

Cpi  =  0 

ACl 

0.0195 

ACp] 

-0.0015 
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